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<~^ . Abstract 

l> 

cn . We investigate the effects of the Kepler rotation of lens binaries on the binary- 

O ' microlensing events towards the Large Magellanic Cloud (LMC) and the Small Mag- 

Q>^ I ellanic Cloud (SMC). It is found that the rotation effects cannot always be neglected 

^ ! when the lens binaries are in the LMC disk or the SMC disk, i.e., when they are 

Q_|. self-lensing. Therefore we suggest that it will be necessary to consider the rotation 

O ' effects in the analyses of the coming binary events if the microlensing events towards 

$—1 ' 

^ [ the halo are self-lensing. As an example, we reexamine the MACHO LMC-9 event, 

in which the slow transverse velocity of the lens binary suggests a microlensing event 

- in the LMC disk. From a simple analysis, it is shown that the lens binary with total 

Vh ' mass ~ IMq rotates by more than ~ 60° during the Einstein radius crossing time, 

c^ ' 

However, the fitting of MACHO LMC-9 with an additional parameter, the rotation 

period, shows that the rotation effects are small, i.e., the projected rotation angle is 
only ~ 5.9°(M/M0)^/^ during the Einstein radius crossing time. This contradiction 
can be settled if the physical parameters, such as the mass and the velocity, are differ- 
ent in this event, the binary is nearly edge-on, or the binary is very eccentric, though 
definite conclusions cannot be drawn from this single event. If the microlensing events 
towards the halo are due to self-lensing, binary-events for which the rotation effects 
are important will increase and stronger constraints on the nature of the lenses will be 
obtained. 
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§1. Introduction 

The analysis of the first 2.1 years of photometry of 8.5 x 10^ stars in the LMC by the 
MACHO Collaboration B^ suggests that the fraction 0.62lo;| of our halo consists of massive 
compact halo objects (MACHOs) of mass O.5]'^o!2^0 ^^ the standard spherical fiat rotation 
halo model. A preliminary analysis of four years of data suggests the existence of at least 
eight additional microlensing events with tjur ~ 90 days in the direction of the LMCB-* 

At present, we do not know what MACHOs are. There have been several identifications 
of MACHOs proposed, such as brown dwarfs, red dwarfs, white dwarfs, neutron stars, pri- 
mordial black holes, and so on.i)'i'0)'i)'i'lil)'i)'i)'§)'0)'§)'i)'0)'§)'S'li)'lii)'§)'i)'§) Any 
objects clustered somewhere between the LMC and the Sun with column density larger than 
25M0pc~^ may also explain the data. ^ They include the following possibilities: LMC-LMC 
self-lensing, the spheroid component, thick disk, a dwarf galaxy, tidal debris, and warping 
and flaring of the galactic disk. ©'B'i) '§)'§)'©'§) (See also Ref. |0D.) 

Such obscurities of the mass and the spatial distribution essentially result from the fact 
that the time scale of an event, which is an important observable, is a degenerate combination 
of the three quantities one would like to know, the mass, the velocity and the position 
of the lensing object. Several methods have been proposed to break these degeneracies, 
for example, launching a parallax satellite into solar orbit, Eil''E3'E3''a'El) observing the 
annual modulation in light magnification induced by the Earth's motion, Ej^'EZP observing 
the deviation of light magnification from a simple point-source model due to the finite- 
source size effect when the impact parameter of the trajectory of the lens is comparable to 
the source size,B'E§'E2P distinguishing the dependence of the lensing rate on the background 
stellar density, cIf and so on. 

A microlensing event due to a binary is one of the best candidates to break the degen- 
eracies. In binary-microlensing events, the light magnification dramatically deviates from 
that of a simple point-source model when the source transverses the caustics, where a point 
source is amplified infinitely. We can obtain information concerning the transverse velocity 
of the source from this deviation, which can be used to distinguish between halo-lensing and 
self-lensing. To this time, two binary-lens microlensing events have been observed, MACHO 
LMC-g© and MACHO 98-SMC-l.Ei'E])'§) Although we cannot say for certain from only 
these events, e3* these events support self-lensing because of slow transverse velocities. In the 
future, the number of binary-lens microlensing events will increase, IllPW and hence these 
events will be important to break the degeneracies in the physical parameters. 

In almost all analyses of binary-lens microlensing events, the rotation of the lens binary 
has been neglected. eHP This is because the period of a lens binary in the halo is much larger 



than the time scale of the amphfication. However, the lensing object of the MACHO LMC-9 
event, for example, is very likely to reside in the LMC disk, not in our halo. Since the 
characteristic transverse velocity of a lens in the LMC disk or the SMC disk is smaller than 
that in our halo, the rotation of the lens binary in the disk may be important. For this reason 
we reconsider the rotation effects on the analyses of the binary-lens microlensing events in 
this paper. As an example we reanalyze the MACHO LMC-9 event taking the rotation into 
account. Note that, in the previous analysis of this event, the transverse velocity is somewhat 
smaller than that expected for a lens in the LMC disk.eS^ If we take the rotation of the lens 
binary into account, the transverse velocity may be larger, since the incident angle of the 
trajectory of the source into the caustics may be smaller, and hence the source may take 
shorter time to move by one stellar radius of the source. This is one of our motivations to 
examine the rotation effects. 

In §2 microlensing by a double point mass is reviewed. In §3 the rotation effects of a 
lens binary are estimated. We suggest the possibility that the rotation effects are important 
when the lens binary resides in the LMC disk or the SMC disk. In §4 we perform the fitting 
of the MACHO LMC-9 event, taking into account the rotation of the lens binary. Section 5 
is devoted to summary and discussion. 

§2. Microlensing by two point masses 

We now briefly review microlensing by a double point massE^'H to introduce our nota- 
tion. We consider a lens binary consisting of two point masses. Mi and M2, whose center 
is at a distance Dd from the observer, and we consider a source at a distance Dos from the 
observer in Fig. |I|. We define the lens plane as the plane which contains the center of mass 
of the lens binary and is perpendicular to the line connecting the observer and the center of 
mass, i.e., the optical axis. We also define the source plane as the plane which contains the 
source and is parallel to the lens plane. The distance between the lens plane and the source 
plane is written as Dig = Dos ~ D^i- We define a coordinate system (^j., ^y) on the lens plane 
and (?7i,., rjy) on the source plane, taking the origin of each coordinate at the intersection 
between each plane and the optical axis. An equation which relates the image position ^ to 
the source position r/ is called a "lens equation" . From Fig. |l], we see that the lens equation 
for lensing by a double point mass is 
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where ^i and ^2 ^^^ the positions of the masses projected onto the lens plane. 0(^) is the 
deflection angle of light due to the lens masses, which is the summation of the deflection 
angle due to each mass. The Einstein radius for the total mass of the binary, M = Mi +M2, 
is defined as 
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where x := Doi/Dos- With the definitions 
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the lens equations (2-1) and (2-2) become dimensionless equations. 
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Note that we are discussing the situation on the lens plane, since we normalize the length 
as in Eq. (p4D. 




observer 

Fig. 1. The geometry of the gravitational lensing by a double point mass lens. 

These equations (2-5) and (2-6) can be used to find all images of the source. tlP If we 
cannot resolve the images by observations, the only observable quantity is the amplification 
of the brightness of the source. The amplification factor I is the inverse of the determinant 
of the Jacobi matrix, 
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Fig. 2. Caustics for a double point lens. The lenses have identical masses and are separated by 
/ = 0.6, 1.4 and 2.2. 

where r j is the image position. For certain values of r, the amplification factor / diverges. A 
set of these points forms curves called "critical curves" . The projection of the critical curves 
onto the source plane with Eqs. (2-5) and (2-6) forms caustics on the source plane. Caustics 
have three kinds of morphology, depending on the separation of the masses. In Fig. 2 these 
three kinds of caustics are shown. The number of images is five in the closed caustics and 
three outside. Around the caustics, large amplification appears, and the light curve has a 
peak. 



§3. Microlensing events by a binary lens 

3.1. Motion oj the source 

We assume that we can neglect the Earth's motion around the Sun, and consider the 
relative motion of the observer, the lens and the source as the motion of the source. Every 
quantity on the source plane is projected onto the lens plane by Eq. (|2-4| ). Hence it is 
convenient to consider the motion of the source as that projected onto the lens plane with 



Eq. (|2-4| ). The trajectory of the source is characterized by its impact parameter h with 
respect to the origin of the lens plane in units of the Einstein radius pe and the angle Q 
between the x-axis on the lens plane and the trajectory. The source closest approaches the 
origin at Ts-, and the transverse velocity of the source in the lens plane is Vj_. The transverse 
velocity of the source pulled back onto the source plane Vt is related to V^ as 



Vt 



X 



The Einstein radius crossing time tE is defined as 



tE-.-- 
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It is convenient to measure the time T with respect to Ie-, t := T /Ie- By normahzing the 
time in this manner, the source moves by a unit length in a unit time. The source closest 
approaches the origin at the normalized time, ts '■= Ts/tE- 

If there are enough data around the peak of the light curve, important information about 
the size of the source can be obtained. eHI* If the radius of the source is R^, it takes a time 

for the source to move by one stellar radius of the source. 

3.2. Rotation of the lens binary 

The lens binary rotates according to Kepler's law. The relative vector U is defined 
by [/ := Ui — U2, where t/i and t/2 denote the positions of the binary masses in the 
orbital plane. We take the origin of the orbital plane at the center of the binary masses as 
M1U1 + M2U2 = 0. Hence, the positions of the binary masses are determined by the relative 
vector U as Ui = 112U and U2 = —^lU, where fii and fi2 are defined in Eq. ( p-4D . In the 
orbital plane, the relative vector U = {U cos (j),U sin (p) at a time T is determined through 
the parameter A as 

T-To = ^(A-esinA), (3-4) 

U = A{l-ecosX), (3-5) 

cos(0 - 0o) = z r, (3-6) 

1 — e cos A 

where 
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is the period, A is the semi-major axis, and e is the eccentricity of the lens binary. Tq and 
00 are integral constants. We can take the coordinates on the orbital plane so that 00 = 0. 
With the definitions 

Tb ^ To U A 

ts ■■=-—, to := — , u= — , a:= — , (3-8) 
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the above equations (3-4), (3-5) and (3-6) become dimensionless: 

t-to = -^(A-esinA), (3-9) 
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Fig. 3. The relation between the orbital plane and the lens plane. The {rx,ry) plane is the lens 
plane, the {ux,Uy) plane is the orbital plane, and a, 13 and 7 are Euler angles. 

In general, the orbital plane does not coincide with the lens plane. Since we consider 
the case Doi ^ A, it is convenient to think about the lens masses projected onto the lens 
plane, i.e., the thin lens approximation. The orientation of the coordinates in the orbital 
plane relative to the coordinates in the lens plane is determined by Euler angles a, (3 and 7. 
In Fig. 3, we show the relation between the orbital plane and the lens plane, where [rxiTy) 
plane is the lens plane, {ux,Uy) plane is the orbital plane, and a, {3 and 7 are Euler angles. 
We can take the coordinates in the lens plane so that a = 0, since a can be absorbed into 
^, the angle between the x-axis in the lens plane and the trajectory of the source. A point 
{u cos 0, u sin 0) in the orbital plane is projected to the point (/ cos v?, / sin ip) in the lens plane, 
where 



tan(/9 



/ = -uy cos^ (3 cos^(0 + 7) + sin^(0 + 7), 
tan(0 + 7) 
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(3-12) 
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3.3. Estimate of the rotation effects 

We can estimate the rotation effects of a lens binary by comparing the period of the 
binary Tb with the Einstein radius crossing time tE- If the ratio ts = Ts/tE, is not much 
larger than unity, we cannot neglect the rotation of the binary. 



For a typical MACHO, the ratio ts is given by 



_Tb_ 27iV^ / A^ _ 2nVs_ p%a^ 

^ \l) l^SOkpcyl \Mq) 1^200 km/syl 1^0.5(1 -0.5) y ' ^ ' 

using Eqs. (|2-3|), (p-2|), ( |3-7| ) and (3-8). This imphes that a typical binary in our halo rotates 
by about 360°/tB ^ 2.1°{M/MqY'^ during the Einstein radius crossing time. This may be 
small enough to neglect the rotation of the binary.lill 

On the other hand, if a lens binary resides in the LMC disk, the parameters are different. 
For a typical LMC lens, assuming that the thickness of the LMC disk is smaller than ~ 500 
pc, the ratio ts is given by 

^~ nj UokpcJ [m^J l^eokm/sAo.ggii-o.gg)^ ' ^^^^' 

which implies that the binary rotates by more than ~ 20°(M/Mq)^/^ during the Einstein 
radius crossing time. Thus, if the lens resides in the LMC disk, the rotation effects may be 
important for the fitting of the light curve. This argument can also be applied to an SMC 
lens, since the transverse velocity is small and x is close to 1. 

Bennett et al. e3) fitted the data of MACHO LMC-9 and obtained fitted parameters as 
in Table I. The radius of the source is estimated as 

i?^ = 1.5 ± 0.2 Rq, (3-16) 

using the theory of stellar evolution. The transverse velocity in the source plane can be 
estimated as 

VV = -^ = 19 ± 6 km/s, (3-17) 

-'* 

with Ti^ in Table I and Eq. ( p-3| ). Comparing this value with the probability distribution of 
the transverse velocity in the LMC disk and the Milky Way halo, a lens in the LMC disk is 
preferred over a halo lens.e^ Assuming that the lens of MACHO LMC-9 resides in the LMC 
disk, the ratio t^ is given by 

*^<^-Kt) (sol^j (mo) (T9k^j(o.99(l-0.99)] ' ^^'^^^ 

where we use the relation V± < Vr, since V± = xVr and x < 1. Since Eq. ( |3-18[ ) indicates 



that the binary rotates by more than ~ 60°{M/MqY^^ during the Einstein radius crossing 

*' For a close binary, i.e., for a <C 1, it seems that the rotation effects can be large from Eq. (3-14). 
However, this is not a correct argument, as discussed in Appendix H. 
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time, we cannot neglect the rotation of the binary for the fitting of the data. Therefore we 
perform the fitting of the MACHO LMC-9 event taking into account the rotation of the lens 
binary in the next section. 

The conclusion of this section is that the rotation effects cannot be always neglected when 
the lens binaries are in the LMC disk or the SMC disk, while the rotation can be neglected 
when the lenses are in the Milky Way halo. 



§4. Fitting of the observed data in the binary-lens events 

4.1. Parameters characterizing the binary-lens microlensing events 

We summarize the parameters necessary to describe the binary-lens microlensing event 
in this section. First, a dual color observation requires two parameters, taking the blending 
into consideration: 

(1) foR '■ fraction of the lensed brightness of the lensed star in the red band. 

(2) foB '■ fraction of the lensed brightness of the lensed star in the blue band. 

When we can neglect the rotation of the lens binary, we need the following parameters: 

(3) tE '■ Einstein radius crossing time. 

(4) Ts : time when the source most closely approaches the origin. 

(5) b : impact parameter of the source in the lens plane in units of the Einstein radius. 

(6) 6 : angle between the x-axis in the lens plane and the trajectory of the source. 

(7) fii : the mass fraction of the first mass. 

(8) / : separation of the lens masses projected onto the lens plane in units of the Einstein 
radius. 

(9) T^ : time for the source to move by one stellar radius of the source. 

If we can neglect the finite size of the source, the last parameter T^ is not necessary. 

When we consider the rotation of the lens binary, we need the following parameters in 
addition to the above parameters: 

(8) a : semi-major axis of the lens binary in units of the Einstein radius. 

(10) tB '■ period of the lens binary in units of the Einstein radius crossing time. 

(11) f3 : Euler angle. 

(12) 7 : Euler angle. 

(13) e : eccentricity of the lens binary. 

(14) To : time when the binary is at the pericenter. 
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where the parameter (8) has been replaced. We need five additional parameters when we 
take the rotation of the lens binary into account. 

If we assume that the orbit of the binary is circular (e = 0), the parameters (13) and 
(14) are not necessary, since To can be absorbed into 7. If we assume that the orbit of the 
binary is face-on, the parameters (11) and (12) are not necessary, since 7 can be absorbed 
into 6. If we assume that the lens binary is face-on and e = 0, the parameters from (11) to 
(14) are not necessary. 

4.2. Fitting of the MACHO LMC-9 event 

We analyze the raw data for the MACHO LMC-9 event. The baseline of the photometry 
that corresponds to no amplification has to be determined by the fitting with it added as 
one more parameter. However, to save time, we determine the baseline by the least squares 
fitting of the data for T < 300 day and 900 day < T, where the amplification is expected 



to be less than (4^ -|- 2)/(4 x y/4P + 4) ~ 1.006. With this basehne we can translate the raw 
data into the data for the amplification. 

We first performed the fitting of the data of MACHO LMC-9 neglecting the rotation of 
the lens binary to check our fitting code. The result of the x^ fit is shown in Table II. lI This 
result does not differ greatly from the result in Table I, which is a crosscheck of our fitting 
code. The structure of the caustics and the trajectory of the source are shown in Fig. 4. The 
corresponding light curves are also shown. 

Next we performed the fitting of the data of MACHO LMC-9 taking account of the 
rotation of the lens binary. The infiuence of the rotation of the lens binary on the light 
curve can be divided into two effects, the change of the caustics shape with time caused by 
the change of / and rotation of the caustics with time in the lens plane, which is treated as 
a curve of the source trajectory. However, since the former effect is hardly separated from 
the ambiguity of the other parameters within the limited accuracy of the observation, we 
consider only the latter effect. Thus, in order to evaluate the rotation effects, it is useful 
to examine the case that the binary is face-on and e = 0. Therefore we assume that the 
binary is face-on and e = as a first step. Then we need ten parameters, (I)-(IO) in §4.1. 
The results of the x^ fitting are shown in Table HI. The structure of the caustics and the 
trajectory of the source are shown in Fig. 5. The corresponding light curves are also shown. 



*' There are generally several local minima in_th£^x£^ fitting using some parameters. Several methods 
have been proposed to find the global minimum.EJ'Ej^'EJ) However, this requires a great deal of effort to 
find the global minimum. Therefore, we restrict the parameter space somewhat from the shape of the light 
curve, and then we pick the parameter set with the smallest x^ from ~ 50 fittings. For example, we can 
restrict the morphology of the caustics to the middle in Fig. 2, since the amplification is sufficiently high 
between the peaks of the light curve. 
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Fig. 4. The structure of the caustics and the trajectory of the source are shown (left). The corre- 
sponding hght curves for red {middle) and blue (right) are also shown. The rotation effects are 
not included. 

Contrary to our expectations, there are few differences between the fitted parameters 
with and without the rotation from Tables II and III. This is because the period of the lens 
binary t^ is quite large. The binary rotates by only ~ 5.9°{M/Mq)^^^ during the Einstein 
radius crossing time. Thus the rotation effects are very small.tJ Note that the result of 
the small rotation effects does not depend on the face-on and e = assumptions. We can 
consider the following three possible reasons for the contradiction of the simple estimate of 
tB in Eq. ( p-15| ) and the fitted value of ts in Table III: 

1. The physical parameters in Eq. (3.15), such as the mass M and the velocity V±, are 
not typical in this event. 

2. The binary of this event is nearly edge-on. 

3. The binary of this event is very eccentric. 

We consider the possible physical parameters that account for the small rotation effects in 
§4.3, we consider the possible inclination in §4.4, and we consider the possible eccentricity 
in §4.5. 

4.3. The mass, the velocity and the distance of the lens 

Assuming that the binary is face-on and that e = 0, we can obtain the mass, the velocity 
and the position of the lens from only the fitted parameters in Table III and the radius of 
the source in Eq. ( p-16 ) as follows. S Note that the probability distribution of the transverse 
velocity is not necessary. 



*■* Of course, the fitted parameters may be at only a local minimum. The period ts in the global 
minimum may be smaller. However, it does not seem that Ib is smaller, since we usually find ts <^ 60 even 
when we start the fitting from is ~ 60. 
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Fig. 5. The structure of the caustics and the trajectory of the source (left). The corresponding 
hght curves for red {middle) and blue (right) are also shown. The rotation effects are included. 



We assume that the binary is face-on and that e = 0. The velocity of the lens V± is 
determined by Eq. (E^) as 



Vi =x- 



R* 



(4-1) 



for given x. Using Eqs. (p-3|) , ( p-2| ) and (|4-1|) , the total mass of the lens M is determined by 



M 



AGD,, \TJ l-x 
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(4-2) 



for given x. Using Eqs. (|2-3|), (|3^), ( |3-7| ), (3-8) and (^4|), the position of the lens x is 
obtained from 



x'^{l — x) 



tl cHeT, 



(4-3) 



The left-hand side of this equation, /(x) := x^(l — x), is plotted in Fig. 6. The function /(x) 
has a maximum value /(2/3) = 4/27 at x = 2/3. The equation /(x) = y has two solutions, 
Xi and X2, for < y < 4/27. The right-hand side of Eq. ( [4 -31 ) is determined from the fitted 
parameters. 

After determining the positions, Xi and X2, with Eq. ( [4 -31) , we can obtain the mass and 
the velocity of the lens for each distance with Eqs. ( ^-11 ) and ( [4 -21 ). Since the right-hand side 
of Eq. (^-31) is 0.0860 from Table III, two solutions of this equation are 



xi = 0.369 and X2 = 0.892. 



The corresponding mass and velocity are 



M = 0.000902 Mq, 14 = 7.30 km/s for x = Xi, 
M = 0.0127 Mq, V± = 17.6 km/s for x = X2. 



(4-4) 



(4-5) 
(4-6) 
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Fig. 6. f{x) = x^(l — x). The horizontal dotted hne is the right-hand side of Eq. (5 • 3) with the 
fitted parameters in Table III. 

If the binary of the MACHO LMC-9 event has such physical parameters, we can explain the 
small rotation effects. However, these parameters seem to be quite strange. The mass seems 
to be too small, and the position favors the halo lens, while the transverse velocity prefers 
the LMC lens. However, we cannot draw definite conclusions from this event alone. 

4.4. Inclination 

In this section we consider a possible inclination to explain the small rotation effects in 
the MACHO LMC-9 event. Without the face-on assumption, the fitted values of ts and a do 
not generally coincide with the real values of ts and a. We obtained the mass, the velocity 
and the position using t^ and a instead of ts and a, respectively, in Eq. ( |4-3| ), where t^ 
and a denote the fitted values of ts and a, respectively. Therefore a certain inclination may 
explain the small rotation effects even if we use typical physical parameters. To investigate 
the inclination effect, we set e = 0. Therefore the additional parameters are the Euler angles 
P and 7, as shown in §4.1. Rigorously, the parameters ts and a for given /3 and 7 have to 
be determined by the fitting. However, we can approximately determine ts and a for given 
f3 and 7 from ts and a as follows. 

Assuming that a is mainly determined by the separation between the projected masses 
at t = ts, the relation between a and a can be estimated by Eq. (3-12) as 



cos^ (3 cos^ 7 + sin^ 7 



(4-7) 



since (p = a.t t = ts when e = 0. To determine the relation between t^ and ts, we assume 
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that the parameter is is mainly determined by the rotation angle projected on the lens 
plane. In other wards, the relation between ts and ts is determined by the condition that 
the projection of the rotation angle Air/tB between t = ts — 1 and t = ts + 1 coincides 
with Ati/Ib- The angle between the relative vector of the binary masses U and the x-axis 
of the orbital plane at t = t^ — 1 is 0i := —l-RJtB- The corresponding angle y^i between the 
relative vector projected on the lens plane and the x-axis of the lens plane is determined by 

Eq. (3-13) as 

tan(0i + 7) 

tanv^i = . (4-8) 

cosp 

Similarly, since the angle between the relative vector of the binary masses TJ and the x-axis 
of the orbital plane at t = ^5 + 1 is 02 = Syr/ts, the corresponding angle ^2 between the 
relative vector projected on the lens plane and the x-axis of the lens plane is determined by 
Eq. (3-13) as 

tan(02 + 7) (. ox 

tanv32 = ^ . (4-9) 

cosp 

The relation between Ib and Ib is determined by the condition \>^2 — '^\\ = ^T^/tB- With 

Eqs. (^4-8|) and ([4 -91) this relation can be obtained as 



27r 

— = arctan 

tB 



=Fcos/3 + Jcos'^ P + cos^7(cos2/5 + tan^ 7) (cos^ /3 tan^ 7 + 1) tan^(47r/tB) 



cos^ 7(cos^ (3 tan 7 + 1) tan(47r/tB) 

(4-10) 

where the minus sign is for < /? < 7r/2 and the plus sign is for 7r/2 < j3 < tt. 

Of course, we cannot determine (3 and 7 by the fitted parameters. However, inversely, 
we can determine the allowed region of (3 and 7 for given x. For example, if we assume that 
the lensing object resides in the LMC disk, i.e., 

0.99 <x<l, (4-11) 

the left-hand side of Eq. (|]3|) is less than 0.009801. Substituting Eqs. (^) and (4-10) into 
Eq. ( [4 -31 ), the allowed region of /? and 7 can be obtained in Fig. |^bJ The probability that the 
inclination and the phase are in the allowed region is 26.0%, assuming random inclination 
and phase. The corresponding mass and velocity can be obtained by Eqs. ( [4- 1]) and ( [4 -21) as 



0.153 Mq < M, 19.6 km/s < ^^ < 19.8 km/s. (4-12) 

In this way, if the binary of this event is nearly edge-on, we can explain the small rotation 

effects with typical physical parameters. However, we cannot draw definite conclusions from 

only this event. 

*' We find that the inequalities x < xi and X2 < x hold for any (3 and 7, since the inequality t\/a^ < 
i%/a'^ holds irrespective of /3 and 7. However, these relations do not hold if we consider the eccentricity. 
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Fig. 7. The allowed region of the inclination (3 and the phase 7 for 0.99 < x < 1. 



4.5. Eccentricity 

In this section we consider a possible eccentricity to explain the small rotation effects in 
the MACHO LMC-9 event. As in the previous section, without the e = assumption, the 
fitted values ts and a do not generally coincide with the real values of ts and a. Since we 
use ts and a instead of ts and a in Eq. ( [4 -31 ) to obtain the physical parameters, a certain 
eccentricity may explain the small rotation effects even if we use typical physical parameters. 
To investigate the eccentricity effect we consider the face-on binary. Therefore, the additional 
parameters are e and Tq, as shown in §4.1. We can approximately determine ts and a for 
given e and Tq from ts and a as follows. 

Assuming that a is mainly determined by the separation between the projected masses 
a,tt = tsi the relation between a and a can be approximated by 

5(1 + e cos 03) .. ^„^ 

a= , 4-13 

1 — e^ 

from Eqs. (3-10) and (3-11), where 03 is determined from Eqs. (3-9) and (3-11) with t = ts- 
To determine the relation between t^ and ts, we assume that the parameter ts is mainly 
determined by the rotation angle projected on the lens plane. In other words, the relation 
between ts and t^ is determined by the condition that the rotation angle between t = t^ — 1 
and t = ts + ^ coincide with An/tB- The angle 04 between the relative vector between the 
masses and the x-axis of the lens plane at t = ^5 — 1 is determined by 



ts 



tn 



-(A4 

2n 

cos A4 



esin A4), 



cos< 



e cos A4 ' 



(4-14) 
(4-15) 
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Fig. 8. The allowed region of the eccentricity e and ts — to (mode t^) for 0.99 < x < 1. 

from Eqs. (3-9) and (3-11). Similarly, the angle 05 between the relative vector between the 
masses and the x-axis of the lens plane at t = t^ + 1 is determined by 



tB 

t5 + 1 - to = --(As - e sin A5) 

ZTT 



COS Ah 



COS(/)5 



1 — e cos A5 
The relation between t^ and t^ is determined by the condition 
Eq. (4-14) from Eq. (4-16), we can determine t^ as 

2n 
T^ 

where X := (A4 + A5)/2 and Y := (A5 — A4)/2 are determined by 

X — e sin X cos Y 

t.Q — tn 



Y — e cos X sin Y, 



(4-16) 
(4-17) 
64 = Air/tB- Subtracting 

(4-18) 



Y — e cos X sin Y 
from Eqs. (4-14) and (4-16), and 

/47r\ e^cos2X + (2 -e^)cos^F - 2ecosXcosy + e^ - 1 



cos 



(4-19) 



(4-20) 



\tB J e^(cos^X + cos^ y) — 2ecosXcosF + 1 — e^ 

from Eqs. (4-15) and (4-17) and the condition 05 — 04 = An/tB- 

We assume that the lensing object resides in the LMC in Eq. ( 4-ll|) . Substituting 
Eqs. ( |4-13|) and ( |4-18| ) into Eq. (|4-3| ), the allowed region of e and ts — to (mode t^) can 
be obtained in Fig. |^. The minimum eccentricity of this allowed region is 0.892. The cor- 
responding mass and velocity are in Eq. (4-12). In this way, if the binary of this event is 
very eccentric, we can explain the small rotation effects with typical physical parameters. 
However we cannot say for certain from only this event. 
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§5. Summary and discussion 

We investigated the rotation effects of lens binaries on the binary-microlensing events 
towards the LMC and SMC. It is found that the rotation effects cannot always be neglected 
when the lens binaries are in the LMC disk or the SMC disk, while the rotation can be 
neglected when the lenses are in our halo. It follows from this that we need to consider 
the rotation effects in the analyses of the coming binary events if the microlensing events 
towards the Magellanic Clouds are self-lensing. 

In the MACHO LMC-9 event the transverse velocity prefers the lens binary in the LMC 
disk.cH' For this reason, we reexamined this event in detail as an example. The simple 
estimate in Eq. (|3-18|) indicates that the binary rotates by more than ~ 6O°(M/M0)^'^ 
during the Einstein radius crossing time. Therefore we cannot neglect the rotation of the 
lens binary in the fitting of MACHO LMC-9. We performed the fitting of MACHO LMC-9 
taking the rotation into account on the assumption that the binary is face-on and that e = 0. 

Contrary to our expectation, the rotation effects are small, i.e., the projected rotation 
angle is only ~ 5.9°(M/M0)^/^ during the Einstein radius crossing time. Note that this 
result does not depend on the face-on and e = assumptions. We can consider three 



possible reasons for this contradiction between the simple estimate of ts in Eq. ( |3-18|) and 
the fitted value of ts in Table III. The first possibility is that the physical parameters in 
Eq. ( p-18| ) are different in this event as in Eqs. ( [4-4|) -(4-6), although these parameters are 
quite peculiar. The second possibility is that the binary of this event is nearly edge-on with 
typical physical parameters 0.99 < x < 1 and Eq.(4-12), as shown in Fig. |^, although the 
probability of (3 and 7 being in the allowed region is only 26% assuming random inclination 
and phase. The third possibility is that the binary of this event is very eccentric with 
typical physical parameters 0.99 < x < 1 and Eq. (4-12), as shown in Fig. |, although the 
minimum eccentricity of the allowed region is 0.892. However, since we cannot determine all 
physical parameters including the inclination and the eccentricity only from the light curve, 
we cannnot draw definite conclusions from only this event. 

The transverse velocity is somewhat smaller than that expected for a lens in the LMC 
disk. This result is the same as that of the analysis of Bennett et al.,E3) even if we take the 
rotation of the lens binary into account. Since the incident angle of the trajectory of the 
source into the caustics is not greatly changed, the transverse velocity remains small. 

The fitted parameters only represent a local minimum. However, the complete analysis is 
difficult, since the number of data points around the peak of the light curve is so small that 
there will be many local minima. E^* The lack of data also allows the possibility that the source 
is a binary, as Bennett et al.E^ claimed. We agree that it is important to determine whether 
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the source is a binary or not with HST. In order not to overlook the important information 
from the caustic crossing events, it is necessary to monitor the events frequently, such as 
with the MOA collaboration!^ and the Alert system. 

Note that the rotation effects are not always important when the lens is in the LMC disk 
or the SMC disk. For example, in the MACHO 98-SMC-l event the transverse velocity is 
larger than ~ 60 km/s.B'cJ'E^ Since t^ ^ 20 ^ 1 from Eq. (|3-15|), the rotation effects 



are not efficient in this event. Our statement is that there is a sufficient probability that the 
rotation effects are important for self-lensing. 

If the microlensing events towards the halo are due to self-lensing, E2P' til) the number of 
binary-events for which the rotation effects are important will increase. We will be able to 
crosscheck whether or not the microlensing events are self-lensing by examining the rotation 
effects, together with the transverse velocity distribution. As the number of binary-lens 
events increases, stronger constraints on the nature of the lenses will be obtained statistically. 
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Appendix A 

Notations 



^ : image position in the lens plane 

77 : source position in the source plane 

^1) ^2 • positions of the masses of the lens binary in the lens plane 

r : image position in the lens plane in units of pe 

z : source position in the source plane in units of pe 

Vi, r2 '■ positions of the masses of the lens binary in the lens plane in units of pe 

U : relative vector between the lens masses in the orbital plane 



t/l, C/2 

u 

Ml, Ms 
M 

/ 

foR 

foB 

Dol 

Dos 

Dls 

X 

PE 

v± 

a 
a 
I 
e 



tE 
t 

Ts 
ts 
Tb 

tB 

Ib 
To 
to 
T. 



position of the masses of the lens binary in the orbital plane 
relative vector between the lens masses in units of pe 
masses of the lens binary 
Ml + M2, total mass of the lens binary 
Mi/M, mass fraction (i = 1,2) 
amplification factor 

fraction of the lensed brightness of the lensed star in the red band 
fraction of the lensed brightness of the lensed star in the blue band 
distance between the observer and the lens plane 
distance between the observer and the source plane 
distance between the lens plane and the source plane 
Doi/Dos 
Einstein radius 

transverse velocity of the source in the lens plane 
transverse velocity of the source in the source plane 
impact parameter of the source in the lens plane in units of pe 
angle between the x-axis and the trajectory of the source 
semimajor axis of the lens binary 
semimajor axis of the lens binary in units of pE 
fitted value of a assuming j3 = and e = 

separation of the lens masses projected on the lens plane in units of pe 
eccentricity of the lens binary 
azimuthal angle of the binary in the orbital plane 

angle between the relative vector projected on the lens plane and the x-axis of 
the lens plane 
Einstein radius crossing time 
T/tE, time in units of tE 

time when the source most closely approaches the origin 
Ts/tE, time when the source most closely approaches the origin in units of tE 
period of the lens binary 
Ts/tE, period of the lens binary in units of tE 
fitted value of ts assuming [3 = and e = 
time when the binary is at the pericenter 
To/tE, time when the binary is at the pericenter in units of tE 
time for the source to move by one stellar radius of the source 
radius of the source 
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a,/?, 7 : Euler angles which determine the orbital plane (see Fig. 3) 

Appendix B 

Rotation Effects of a Close Binary 



In this appendix, we show that the rotation effects of a typical MACHO binary are small 
even when the binary is the close one. For a close binary (i.e. a -C 1), it seems that the 
effects of the rotation can be large from Eq. (3-14). For example, ts is about 4.2 when a is 
0.1, cJ which implies that the binary rotates by about 85° during the Einstein radius crossing 
time. However, we have to note that the region where the rotation effects are important is 
not within the Einstein radius but within the caustics. This is because far from the caustics, 
the light curve due to two point masses is almost the same as that due to a point mass lens. 
Therefore it is not reasonable to compare the binary period with the Einstein radius crossing 
time in the estimate of the rotation effects in Eq. (3-14). Instead, we have to compare the 
binary period with the caustics crossing time. 

To treat the problem analytically, we assume that the mass ratio is unity: /ii = /i2 = 1/2. 
When the semimajor axis is small (i.e. a — > 0), the size of the caustics Za^n satisfiesEB) 



a2 



(B-1^ 



Thus, when the semimajor axis a is small, the ratio of the binary period to the caustics 
crossing time, t^*^, is approximately 



Tb 2 



j.cau ^ J. 



B ~ z — :~ ~ ~}B 

1/4 / , . X -1/4 / ,. X / ,^ X X 1/4 



;a\-yU Dos y ( M\~''" ( V^ \( xil-x) 
2 X 134 ' ^ / 06 \ / \ / ± \ / V ; 



,1/ VSOkpc/ \MqJ V200 km/s/ V0-5(1-0.5)/ 
oc a-^/2 (B-2) 

for a typical MACHO. The exponent with respect to the semimajor axis a changes from 3/2 
to —1/2. Therefore the rotation effects are small even when the binary is the close one, as 
long as the binary is in the Milky Way halo. 
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Table I. Fitted parameters for MACHO LMC-9 from Bennett et al. (1996). The rotation effects 
are not included. 



parameters 


fitted values 


foR 


0.259 ±0.002 


foB 


0.174 ±0.001 


Ie [days] 


71.7±0.1 


Ts [days] 


603.04 ±0.02 


b 


-0.055 ±0.001 


9 [radian] 


0.086 ±0.001 


/ii 


0.620 ±0.002 


/ 


1.6545 ±0.0008 


T. [days] 


0.65 ±0.18 


X^(for 848 degrees) 


1489 


reduced x^ 


1.76 



Table II. Fitted parameters for MACHO LMC-9 from our fitting code. The rotation effects are 
not included. 



parameters 


fitted values 


foR 


0.255 


foB 


0.176 


Ie [days] 


71.4 


Ts [days] 


603.1 


h 


-0.0538 


9 [radian] 


0.0843 


/ii 


0.617 


a 


1.66 


T. [days] 


0.611 


X^(for 871 - 9 degrees) 


1396 


reduced x^ 


1.62 
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Table III. Fitted parameters for MACHO LMC-9 from our fitting code. The rotation effects are 
included. 



parameters 


fitted values 


foR 


0.246 


foB 


0.172 


Ie [days] 


69.4 


Ts [days] 


601.4 


h 


-0.0319 


9 [radian] 


0.144 


/Xl 


0.560 


a 


1.69 


T^ [days] 


0.611 


tB 


60.6 


X^ (for 871- 10 degrees) 


1372 


reduced x^ 


1.59 
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